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ABSTRACT 

We analyse a controlled A-body + smoothed particle hydrodynamics simulation 
of a growing disc galaxy within a non-growing, live dark halo. The disc is continuously 
fed with gas and star particles on near-circular orbits and develops a bar comparable 
in size to the one of the Milky Way (MW). We extract line of sight velocity v\os distri¬ 
butions from the model and compare it to data recently obtained from the APOGEE 
survey which show distinct high velocity features around i;ios ~ 200kms“^. With an 
APOGEE like selection function, but without any scaling nor adjustment, we find uios 
distributions very similar to those in APOGEE. The stars that make up the high wios 
features at positive longitudes I are preferentially young bar stars (age r < 2 — 3 Gyr) 
which move away from us along the rear side of the bar. At negative I, we find the 
corresponding low ujos feature from stars moving towards us. At ( > 10 degrees the 
highest ?;ios stars are a mixture of bar and background disc stars which complicates 
the interpretation of observations. The main peak in uios is dominated by fore- and 
background stars. At a given time, ^ 40 — 50 per cent of high z;ios stars occupy Xi 
like orbits, but a significant fraction are on more complex orbits. The observed feature 
is likely due to a population of dynamically cool, young stars formed from gas just 
outside the bar and subsequently captured by the growing bar. The high ujos features 
disappear at high latitudes |6| > 2 degrees which explains the non-detection of such 
features in other surveys. 

Key words: Galaxy: bulge; methods: numerical - Galaxy: evolution - Galaxy: kine¬ 
matics and dynamics - Galaxy: structure; 


1 INTRODUCTION 


The central region of our Galaxy is dominated by a bar 
(Blitz & Spergel 1991| Binney et al. 19911, which consists 
of an X-shaped box/peanut bulge at i? < 2kpc (Wegg & 


Gerhard]|2013 1 and a vertically thin part, the long bar, ex¬ 
tending to 7? ~ 4 — 5kpc (Wegg et al.||2015l. Recently, the 


Apache Point Observatory Galactic Evolution Experiment 
(APOGEE) survey ( Allende Prieto et al. [2008 1 has revealed 
that the line of sight velocity uios distributions of ~ 4 700 
K/M giant stars towards the Galactic centre at Galactic co¬ 
ordinates Z = 4 — 14 degrees and |6| ^ 2 degrees ex hibit dis- 
tinct shoulders at high velocities vios ~ 200kms“^ (Nidever 
et al.|[2012 I . This result was followed up by other surveys. 


which yielded differing results: On the one hand, |Babusi-| 
aux et al. (20141 analysed the uios distribution of ~ 400 bar 
red clump stars at I = (—6) — 10 degrees and |&| ^ 1 de¬ 
grees and found hints of a distinct high uios component at 
wios > 200kms“^ at positive Z and additionally of a distinct 
low vios component at vios < —200kms“^ at negative Z. On 
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the other hand, Zoccali et al. (20141 failed to detect dis¬ 
tinct high or low wioa components in a survey of ~ 5 000 red 
clump stars in the central Galaxy. However, they were prob¬ 
ing different sight-lines, most of them at |6| >2 degrees and 
thus further away from the plane than the other authors, or 
pointing more closely towards the Galactic centre at Z = 0. 

Bars are rotating triaxial structures which are expected 
to be supported mostly by stars on prograde orbits parented 
by the xi family of closed long-axis orbits ( [Contopoulos 
Papayannopoulos| 1980|). These orbits have proven essential 


in explaining the (Z,u) diagram of gas towards the Galactic 
centre ( [Binney et al.|[T99ip and are thus expected to be an 
important ingredient for explaining any observation of the 


kinematics of bar stars. Molloy et al. (20151 have recently 


shown that they could recover a distinct high uios peak in an 
A-body simulation of a barred galaxy if they analysed only 
2:1 resonant orbits elongated along the long bar axis, i.e. xi 
orbits. They did not recover the feature when taking into 


account all stars. Previously, Nidever et al. (20121 and Li et 


al. (20141 had failed to recover the observational feature in 


A-body simulations of Milky Way (MW) like galaxies. This 
led Li et al. (20141 to conclude that the high uios peak is 
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either a spurious feature due to low numbers of stars or not 
connected to stars on bar orbits. 

Here we analyse a controlled simulation of a growing 
MW-like disc galaxy with respect to nioa distributions to¬ 
wards the barred galactic centre. The model galaxy features 
a live dark halo and is continuously fed with stellar and gas 
particles on near circular orbits. This type of simulation was 
introduced by Sellwood & Carlberg (19841 for a 2D disc but 
has scarcely been used since (e.g. Berrier fc Sellwood||2015 


for a specific problem). The disc develops a bar early on 
and becomes comparable in structure to the MW at later 
stages, which allows us to create uios distributions and com¬ 
pare them to the APOGEE data. 

With the help of this simulation, we explore the origin 
of the high uios structures and ht the detections and non¬ 
detections into a consistent picture. We find that these fea¬ 
tures are associated with a kinematically cool population of 
stars running alongside the bar and likely recently captured 
from the Galactic disc as a consequence of bar growth. 

Our paper is organised as follows: Basic information 
about the data and simulation are found in Section 2 (survey 
data, selection function) and Section 3 (simulation setup). 
The general evolution of disc and bar in our simulation is 
discussed in Section 4. Detailed kinematics are analysed in 
Section 5 with focus on the high uios peak. In Section 6 we 
apply the APOGEE selection function and directly compare 
the model to the data. We also explore some variation of as¬ 
sumed parameters, e.g. latitude h and bar angle 4>- In Section 
7, we use the successfully identified high wios peak for a de¬ 
tailed study of the contributing orbits. Finally, we conclude 
in Section 8 . 


OBSERVATIONS AND SELECTION 
EUNCTION 


For kinematic data from the central bar/bulge region of 
the MW we use data from the Apache Point Observatory 


Galactic Evolution Experiment (APOGEE, Majewski et al. 


2015), an infrared, high-resolution survey, which is part of 
the Sloan Digital Sky Survey III (SDSS [Eisenstein et alT 


2011). The fibre-spectrograph allows for 300 spectra per 
plate covering a spectral range of 1.51 — 1.7/im in the H 
band. As an infrared survey, APOGEE has a unique capa¬ 
bility to obtain information from stars in the high-extinction 
regions within the Galactic plane and towards the central 
region of the MW. We use the APOGEE data publicly 


available via the CAS for SDSS DRIO (Ahn et al. 2014). 


Throughout this work, we correct all observed stellar veloc¬ 
ities from APOGEE to the Galactic standard of rest, us¬ 
ing {U,V,W)sun = (14.0,250.0,6.0), in concordance with 
Schonrich (2012) and [McMillan (2011). 

Another advantage of the survey is the benign selec¬ 


tion function described in Zasowski et al. (2013). Stars 


in the target magnitude range with dereddened colours 
J — Ks > 0.5 mag are targeted evenly. In the case of the cen¬ 
tral Galaxy fields used within this work, stars were targeted 
for 6 < 77 < 11. To model the selection function we employ 
the population synthesis machinery from the [Schonrich fc| 
Binney (2009) model. We apply a standard Salpeter IMF 


with an exponent of 7 = —2.35, neglecting binarity of stars. 
The population synthesis uses a dense grid of B.A.S.T.I. 


selection function 
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Figure 1. The Selection function in distance s and age r at a fixed 
metallicity of Z = 0.0217. Top panel: Contour plot as a function 
of s and t. Bottom panel: The age dependence at a typical bar 
distance of s = 7.5 kpc. Note the strong bias towards stars with 
ages < 2 Gyr at larger distances. 


stellar models ( [Pietrinferni et al.|2004| |2009| for their colour 
calibrations see Bedin et al. 2005). At each distance s and 
age T we count the number of stars per mass of the ini¬ 
tial population that would end up within the colour and 
magnitude selection of APOGEE. The result is shown in 
Figure for a fixed metallicity oi Z — 0.0217. The use of a 
single, slightly super-solar metallicity is a simplihcation, as 
the inner MW has a relatively wide metallicity distribution 
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around this value ( Hill et al.|2011 Bensby et al.|2013 1. How¬ 
ever, all relevant metallicities show a similar preference for 
the detection of young stars. As the detailed age-metallicity 
distributions for the relevant sightlines are not well known 
and the detailed connection between model galaxy and MW 
in terms of structure and formation history is unclear, we 
choose to keep the model simple and use one typical metal- 
licity. The band of high selection probabilities at very small 
distances derives from main-sequence dwarf stars. The stel¬ 
lar models are not fully reliable in this region for such cool 
dwarfs. However, this foreground is strongly suppressed by 
the geometric factor (s^) of the pencil beam survey. The 
number of stars currently in the giant stages of interest to us 
declines strongly (by roughly one order of magnitude) with 
age, implying a strong bias in favour of young stars within 
the APOGEE survey. We apply this selection function to 
our A^-body simulation by weighting each particle with its 
respective selection probability. APOGEE uses dereddened 
magnitudes for targeting, however there are residual effects 
from reddening uncertainties, which we neglect here. 


3 SIMULATIONS 


In this paper we present a controlled collision- 
less+hydrodynamic simulation of a growing galactic 
disc embedded in a non-growing dark matter halo. The 
simulations were all carried out with the Tree Smoothed 
Particle Hydrodynamics (TreeSPH) code GADGET-3, last 
described in Springel (20051. The applied gravitational 


softening lengths are edisc = 30 pc for stellar and gaseous 
disc particles and edm = 100 pc for halo particles. We use 
an opening angle for the tree code of 0 = 0.5. An adaptive 
time-stepping scheme which assigns time steps in bins j ac¬ 
cording to Ati = t/2 ^ < ^J2r]eil |ai| to individual particles 
is applied. Here rj — 0.02 is an accuracy parameter and 
is the gravitational acceleration at the position of particle 
i with softening Ci. The base timestep r is appropriately 
chosen by the code and a maximum of 10 Myr is allowed 
for timesteps. The hydrodynamical timestep is based on a 
Courant-like condition oc nhijca, where k = 0.15 is 

the Gourant parameter, hi is the smoothing length and Cs is 
the soundspeed. We apply an isothermal equation of state 
P — pCg with Cs = 20kms“^ and use Asph = 48 neighbours 
for the smoothing kernel. For further code details we refer 


to Springel (20051. 


3.1 Initial Conditions 

To generate initial conditions (ICs) for our numerical ex¬ 
periments we use the publicly available GallG code 
, which produces near-equilibrium 
collisionless systems with given density 
distributions using an iterative approach. 

The initial system consists of a dark halo with a mass 
of Mdm = 0.995 X 10^^ Mq represented by Adm = 1 x 10® 
particles and an embedded thin galactic disc with a mass 
of Mdiac.i = 0.5 X lO’^® Mq represented by Adisc.i = 5 x 
10® particles, so that the disc particle mass is rridisc = 1 x 
10^ Mq. 


& Springel 2014 


multi-component 


(Yurin 
ICs of 


The dark halo has a Hernquist profile with 


PDM(r) 


Mdm a 
271 r (r -I- a)® 


( 1 ) 


The inner profile is adjusted so that it is similar to an 
NEW profile with concentration c = 9 and virial radius 
R 20 Q = 162.6 kpc. This leads to a scale radius a = 30.25 kpc. 
The halo initially has a spherical density profile and radi¬ 
ally isotropic kinematics, i.e. equal velocity dispersion in the 
principal directions, ar = 04 , = ere and consequently a value 
of the anisotropy parameter of /3 = 1 — (erg -1- cr^) / (2 ct^) = 
0 . 

The disc is set up with a mass profile 


Pdisc.i (A, z') 


Mdisc.i 1,2 f Z 

-Gysech — 

iTTh^hii \hz 



( 2 ) 


with an exponential scalelength of hu = 1.5 kpc and an 
isothermal vertical profile with a radially constant scale- 
height of hz = 107 pc. The vertical velocity dispersion 
thus declines with radius. We assume {o-zjau)^ = 0.5, so 
that Toomre’s Q shows a minimum value of Qmin = 1.15 at 


R ~ hji. 

For further information regarding the ICs we refer to 
Aumer et al. (in prep.). We note that here the halo is rather 
poorly resolved with Adm = 1 x 10® and a rather short soft¬ 
ening length (thaio = 100 pc) is applied. The low Adm is 
due to limited computational resources and a preference to 


use them to resolve the disc. As expected (e.g. Dubinski et 
al.||2009 1, purely collisionless (no SPH) tests with a higher 


Adm = 5 X 10® show reduced disc heating and less bar slow¬ 
down. Yet the disc heating for the simulation studied here is 
comparable to the MW, with similar velocity dispersions for 
10 Gyr old solar neighbourhood stars, since heating by mas¬ 
sive halo particles compensates the lack of heating by Giant 
Molecular Clouds or other sources such as satellite galaxies. 
As far as bar slowdown due to angular momentum transfer 
to the halo is concerned, the radius and rotation frequency 
of bars differ by less than ~ 20 per cent at late times (6-10 
Gyr), the time relevant for this work. 


3.2 Feeding the disc 


To model the continuous growth of galactic discs via star 
formation and gas accretion we, over a timespan of 10 Gyr, 
add new particles with mdisc = 1 x 10'* Mq to the existing 
disc every 5 Myr. We assume a mass growth rate history 
given by 

M(t) = Mo X exp r-—(3) 

\ ^feed / 


with an exponential decay timescale tfeed = 8 Gyr, which 
is motivated by the findings of [Aumer &; Binne^ ( |2009| ), 
and a normalisation Mo = 8 Mq /yr, so that a total of ~ 
4.5 X 10^® Mq are fed over 10 Gyr. The added disc mass 
is radially distributed as p{R) oc exp {—R/hii{t)), with an 
exponential scale-length that increases with time as 


hR{t) =1.5-1- 


(h_B,final — 1.5) / t 


kpc. (4) 


^10 V IGyr^ 

This increase from hn = 1.5 kpc to a scale-length at 10 Gyr 
hRfinal = 4.3 kpc is motivated by the findings of jBovy et al.j 
(20121. The particles are randomly scattered over azimuth 
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bar radial extent 


bar rotational frequency 


circular velocity curve 





Figure 2. The evolution of bar parameters with time: Left-. The bar radial extent Rban defined as the radius where the logarithm of the 
m = 2 Fourier amplitude ln(A 2 ) drops below -1.5. Middle: The rotational frequency of the bar flbar- Right: The azimuthally averaged 
circular velocity curve in the plane of the model galaxy at 7.3 Gyr (solid) and 9.3 Gyr (dashed). The contributions from the 

dark halo are shown in red and the contributions from the baryons are shown in blue. 


(j> and placed at 2 = 0. The coordinate system is regularly 
updated to be centred on the centre of mass of the system. 

The particles are assigned near circular orbits. We de¬ 
termine the new particles’ rotational velocities as = 
\/aR,iRi for a particle i at cylindrical radius Ri and the com¬ 
ponent of the gravitational acceleration pointing towards the 
centre, a_R,i. For stellar particles, we add random velocity 
components in all three directions rj), R and z, drawn from 
Gaussian distributions with a = 6kms“^. 

After a disc goes bar unstable, there are no stable cir¬ 
cular orbits in the bar region. So we impose an inner cutoff 
radius Rent, within which no particles are introduced. We 
use the radius where the logarithm of the m = 2 Fourier 
Amplitude ln(A 2 ) drops below —1.5 as an inner cutoff ra¬ 
dius. 


3.3 The gas component 

In disc galaxies, bars slow down and grow in size due to the 


transfer of angular momentum to the dark halo (Debattista 
fc SellwoodpOOO I. This loss is to some extent balanced by 
the transfer of angular momentum from the gas to the bar 
I Berentzen et al.|2007 1. Therefore, it is desirable to include a 
gas component in barred disc galaxy models. We do not wish 
to introduce a full model containing the accretion of gas onto 
the galaxy, its cooling and the subsequent star formation and 
feedback processes (see e.g. Roskar et al.|20l'2 for an isolated 
galaxy model), as uncertainties regarding the modelling of 


the physical processes are high (see e.g. Fraternali & Binney 
2008 Nelson et al.||2013 l. 

To keep the model simple, we add a constant fraction 
of gas to our galactic disc and model it with an isothermal 
equation of state P = pej with Cs = 20kms“'^. We simulate 
the growth of the disc by continuously adding stellar and 
gas particles. The initial gas disc is created by turning 5 per 
cent randomly chosen stellar particles into SPH particles. 
This mild disturbance of the ICs is unimportant compared 
to the rapid onset of disc instability in the growing galaxy. 
During disc growth we aim at a gas mass fraction in the 
disc of 10 per cent as observed in modern day MW like 


disc galaxies. We achieve this by adding 10 per cent of new 
particles as gas when the mass fraction of gas in the galaxy 
is above 10 per cent and 20 per cent of new particles as gas 
whenever it is lower. 

Computationally, high gas densities in the central 
galaxy are very expensive. Therefore we include a num¬ 
ber of measures to limit the central density: i) Prior to the 
bar forming and the inner cutoff on particle insertion com¬ 
ing into play, we create no gas particles at P < 1 kpc. ii) 
We include a model for star formation in the very central 
galaxy. We turn SPH particles with hydrogen number den¬ 
sities n > nth = 10 cm“® and specific angular momentum 
3 < jth = lOOkms”^ kpc into star particles with a probabil¬ 
ity p = rjtdyn /where t; = 0.1 is an efficiency parameter, 
idyn = 1 / 4:7:0p is the local dynamical timescale and Ati 
is the timestep of the particle. To mimic the effect of a cen¬ 
tral galactic outflow (Bland-Hawthorn & Cohen 20031 we 
only turn one out of five particles into a star and remove 
the other 4 out of five from the simulation. Note that only 
gas particles in the inner bar region are affected by these 
prescriptions. 


4 STRUCTURE OF THE BARRED DISC 


In this section we discuss the structure of our model of a 
barred disc galaxy. Note that we are comparing a specific - 
and in no way fitted - model to a specific galaxy, the Milky 
Way. Though our model should not provide a perfect match 
to our Galaxy, we expect the model to qualitatively share 
important characteristics. Note that due to computational 
limitations we currently have only one model of this type 
including hydrodynamics available. 

In Figure we present the evolution of bar parameters 
in our simulation. We determine a bar radial extent Pbar by 
first calculating the m = 2 Fourier amplitude profile 


N{R) 


■<»(«) -AmY. 




(5) 


NiR) ^ 

where N (R) is the number of particles in the radial bin cen- 
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Figure 3. Stellar surface density S projections of our galaxy model at 7.3 Gyr. Top Left: A Face-On view of the bar region. Top Right: 
An Edge-On view of the bar region. Only stars within 2 kpc of the x — z plane are considered. Bottom Right: A Face-On view of the 
stellar disc. Bottom Right: A Face-On view of the gas disc. The gas surface density was enhanced by a factor 10. 


tred on R and we only take into account stellar particles. 
Then we define /?bar as the radius where \nA 2 {R) drops be¬ 
low -1.5 or as zero if it never reaches the value in the central 
area. 


In the left panel of Figurej^we see that the bar forms in 
the first Gyr, but its radius and strength are varying signif¬ 
icantly. From t ~ 1 Gyr on the bar grows almost linearly in 
time to reach i?bar ~ 4 kpc at t ~ 9 Gyr. Our definition for 
7?bar allows no exact comparison to the values cited for the 
length of the MW long bar. Wegg et al. (20151 determine 
a length of 5 kpc and give 4 kpc as a lower limit, not in 
disagreement with our model. 


In the middle panel of Figure we present the evo¬ 
lution of the bar rotational frequency Ilbar determined by 
following the angle of the bar major axis from snapshot 


to snapshot. The bar initially rotates with a frequency of 
60 — 80kms“^ kpc“^, but as the loss of angular momentum 
to the halo is not balanced by the transfer of angular mo¬ 
mentum from the inflowing gas component, flbar steadily 
declines to reach ~ 25kms”^ kpc“^ at the end of the sim¬ 
ulation. For the MW bar this quantity is little constrained. 
While our value is lower than most estimates for the MW 


| Dehnen|2000 Antoja et al.|20i4 Sormani et al.|2015|), the 

value is in concordance with a recent estimate by |Portail et| 
al. (2015a I (see their Section 8.3 for a discussion of deduced 
MW bar pattern speeds). Our model is thus at the lower 
end of values inferred for the MW. 


In Figure we present surface density projections of 
the whole galaxy model and the bar region at a simula¬ 
tion time t = 7.3 Gyr. The lower panels show face-on views 
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Figure 4. Line of sight velocity njog distributions Al(uiog)/7Vniax of stars at < 4kpc as seen from R = 8.3 kpc at a bar angle of (/> = 20 
degrees for eight different viewing longitudes 1 = 14 — (—10) degrees. We consider stars at latitudes |fe| < 1 degrees and all panels cover 
2 degrees width in 1. We show the distributions for all stars (blue) and stars younger than 2 Gyr (red). The uiog bins used here are all 
8 km s"'^ wide. 


of the stellar and the gaseous components. Apparently, the 
galaxy is dominated by the bar within i? = 4 kpc and shows 
flocculent spiral structure at outer radii. The upper panels 
give close ups of the bar region. The face-on view shows 
two distinct radial regions in terms of minor axis extent. At 
R < 2 kpc, the bar shows a mildly elongated structure, at 
R — 2 — 4 kpc a thinner, long bar structure is visible. The 
edge-on view reveals, that the region R <2 kpc is dominated 
by an X-shaped structure with the tips of this structure at 
{x, z) ~ (2,1.3) kpc, very similar to the structure of the MW 
bulge/bar region inferred by Wegg & Gerhard (20131. This 
view also shows that the outer bar is significantly thinner. 
Around the tips of the bar we measure a vertical exponential 
scale-height of 180 pc, similar to the long bar scale-height in¬ 
ferred by Wegg et al. (20151, but without a super-thin 45 pc 
component as found by these authors. Such a component is 
likely too thin to be resolved in our simulation, considering 
our softening length of 30 pc. 

At a solar neighbourhood like radius R = 8 kpc, we find 
that our model at 7.3 Gyr has a stellar exponential scale- 
length of 3.5 kpc, which grows mildly towards 10 Gyr. The 
disc is thus likely more extended than that of the MW, for 
which estimates usually range at 2-3 kpc. The model has 
a total baryonic mass of 4.2 x 10^® Mq at 7.3 Gyr, it has 
grown from a mass of 0.5 x 10'^° Mq at t = 0. 

In the right panel of Figurej^we plot the circular veloc¬ 
ity curve VciTc{R) in the plane of the model galaxy at 7.3 and 
9.3 Gyr. The contributions from dark matter and baryons 
are shown in red and blue. We measure Ucirc = '/orR by 
calculating the component of the gravitational acceleration 
pointing towards the centre an on an {R, R) grid at 2 = 0 


and at each R average over the equally distributed values 
of <p. We note that the circular velocity curve in the centre 
of the MW is not well constrained and that it is uncertain 
how the observations of gas velocities in the MW bar are 
connected to our definition of Vc.iic{R)- There is only mild 
evolution from 7.3 to 9.3 Gyr visible in Ucirc(7?). 

The circular velocity curves show features at the bar 
radius, deriving from the baryonic matter distribution. The 
reason for these features is that the bar redistributes matter 
towards the centre, whereas by construction of the model 
new baryons are added outside the bar region, which gives 
rise to the baryonic dip at 7? ~ 3 kpc. In our model, the bary¬ 
onic contribution to the circular velocity is greater than the 
dark matter equivalent only within R ~ 2.5 kpc. It is un¬ 
clear how this relates to the MW, but our model agrees 
with two important constraints: A) At i? = 8 kpc and 7.3 
Gyr we recover Udrc = 236kms“^, in good agreement with 
what |Schonridi| ( |2012| finds from local kinematics. B) The 
total dark matter mass within a sphere of r = 8 kpc is 
6.3 X 10^° Mq in good agreement with the dark matter mass 
within a sphere of solar galactic radius that |Piffl et al.| ( f2014[ ) 
have deduced. 


At the evolution stages relevant for this paper, the 
model has a slow bar with Rba.i ~ 0.5i7cR, where Rcr is 
the bar corotation radius. This is likely connected to the fact 
that our galaxy model is sub-maximum, i.e. the dark matter 
contribution to the circular velocity curve is significant in the 
bar/disc region, as fast bars with 7?bar = 0.7 — 1.07?cr are 
thought to require maximum discs with baryon-dominated 
circular velocity curves (Debattista fc Sellwood|[2000 1. 
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Figure 5. Line of sight velocity distributions -/V(uios)/^max of young stars with ages < 2 Gyr at R < 4kpc as seen from R = 8.3 kpc 
at a bar angle of </> = 20 degrees for eight different viewing longitudes i = 14 — (—10) degrees. We consider stars at latitudes |fe| < 1 
degrees and all panels cover 2 degrees width in 1. We show the distributions for stars at three different times: 7.3 (blue), 8.3 (red) and 
9.3 Gyr (green). The uioa bins used here are all 8kms“^ wide. 


5 t;ios DISTRIBUTIONS 

In the previous section we have shown that the bar emerging 
in our galaxy model has comparable structural parameters 
to the MW bar. We now extract line of sight velocity uios dis¬ 
tributions from our model and qualitatively compare them 
to the results of Nidever et al. (20121. In this Section, we 
limit ourselves to the general understanding of the observed 
phenomena. We will attempt a quantitative comparison to 
observed data in the following Section. 


5.1 uios distributions for the bar region 

As we expect stars moving along the major axis of the bar to 
be the source of the stars that comprise the high uios peak, 
we focus here on the central region of the model, dehned 
by galactocentric R < 4 kpc. The full model including the 
selection function will be discussed in the next Section. 

To create uios distributions we place ourselves at a solar 
radius ot R — 8.3 kpc ( |McMillan|201ip and orient the bar to 
an angle of <!> between its major axis and the Sun-Galactic 
centre sightline. Values found for (j) in the literature range 
between 15 and 40 degrees. For most of this paper we choose 
<j) = 20 degrees and show in the next Section that our con¬ 
clusions are unchanged for the full range of cited values. We 
usually define sightlines which are 2 degrees wide in I and 
b and centred on the given values. To curb Poisson noise 
we use the tt (m = 2) symmetry of the model and use two 
sightlines per snapshots, 180 degrees apart in azimuth. For 
each plot we stack 50 snapshots, separated in time by 1 Myr 


and rotate the bar to achieve the same value of (f) for each 
snapshot. 

In Figure]^ we show vios histograms at latitudes |6| < 1 
degrees at a simulation time of 7.3 Gyr for eight different 
longitudes I between 14 and —10 degrees. The figure shows 
histograms for all bar stars with R < 4 kpc in blue and for 
the subset of young star particles with ages < 2 Gyr in red. 
The uios distributions for all stars (blue lines) do not reveal 
signihcant features at high velocities. This is consistent with 


recent examinations of simulations in Nidever et al. 120121, 
I Li et ah] ( |2014[ ) and [Molloy et al.| ( |2015[ ), who also detected 
no features when all stars are taken into account. 


Nidever et al. (20121 reported peaks of high velocity 
stars with uios ~ 200 — 250 km s”'^ of varying extent for lon¬ 
gitudes I between 14 and 4 degrees and [Babusiaux et ah] 
(20141 reported such peaks for I — 10 and 6 degrees and a 
peak at uios ~ —225kms“^ at I = —6 degrees. Gloser in¬ 
spection of the blue histograms at I = ±6 and ±10 degrees 
reveals an excess of stars at the expected uios values. More¬ 
over, the position of the main peak in the Uios distribution 
decreases with decreasing I as observed by [Nidever et al.| 
(20121. These main peaks of the distributions are, however. 


wider in the model than in the observation. 


The fact that the APOGEE selection function (Eig- 
ure[^ favours young stars motivates to examine only stars 
with ages < 2 Gyr shown as red curves. Note that stars 
formed from low angular momentum central gas do not 
play a role for this analysis and are excluded here. The 
main uioa peaks are much narrower for young stars than 
for all stars, indicating that these populations are cooler. 
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Figure 6. Line of sight velocity njoa distributions Af(uios)/^max of stars at < 4kpc as seen from R = 8.3 kpc at a bar angle of iji = 20 
degrees at longitudes Z = 6 it 1 degrees and latitudes |fe| < 1 degrees, uiog bins used here are all 8kms“^ wide. Panels A-D: Populations 
of stars at four fixed age bins, < 1 Gyr, 1—2 Gyr, 2-4 Gyr and > 4 Gyr at three different simulation times: 7.3 (blue), 8.3 (red) and 9.3 
Gyr (green). Panel E: Populations of stars at the four age bins < 1 Gyr (black), 1-2 Gyr (blue), 2-4 Gyr (green) and > 4 Gyr (red) at 
simulation time 7.3 Gyr. Panel F: The population of stars born between 6.3 and 7.3 Gyr in the simulation as viewed at 7.3 (blue), 8.3 
(red) and 9.3 Gyr (green). 


The red histograms clearly show distinct components at 
Pios = 200 — 300kms“^ for Z = 3,6 and 10 degrees and 
at Hies = —300 — (—200) kms“^ for I = —3, —6 and —10 de¬ 
grees, in agreement with observations. Interestingly, we do 
not recover a separate component at Z = 14 degrees, where 


Nidever et al. (2012 I also claimed to find a distinct high vios 


peak, a finding which we will return to in the next Section. 


5.2 Evolution with time and age 

As the model young populations show the observed signa¬ 
ture, it is interesting to ask whether this is specific to a 
unique time or stage of bar evolution or a generic feature of 
young stars. In Figure we examine the evolution of young 
stellar populations with time plotting the populations which 
at respective times are less than 2 Gyr old. We consider sim¬ 
ulation times 7.3, 8.3 and 9.3 Gyr. We notice several points: 
i) There is very little change in the general shape of the 
histograms and in the position and width of the main uios 
peaks at all longitudes Z. ii) The extreme values of uios in¬ 
crease with time by around 20kms“^ for all sightlines which 
do not cross the galactic centre, iii) The high |wios| peaks 
identified in Figure [^become slightly more distinct in time 
and also move to higher velocities, iv) For 9.3 Gyr additional 
features/peaks become apparent at the opposite side of the 


main peaks, e.g. at —150 < uios/kms ^ < —50 at Z = 6 
degrees or at 0 < uios/kms“^ < 100 at Z = — 6 degrees. 

The high uios features at positive longitudes and low 
Wios features at negative longitudes are thus stable over time 
and not connected to any temporary phenomenon. Due to 
the lack of opposite peaks in the APOGEE data we will 
focus our analysis on the simulation time 7.3 Gyr. 

So far we have shown that young stars show high |uios| 
peaks, whereas the whole population of bar stars does not. 
It is thus interesting to study in more detail how the uios 
distribution evolves with age. We attempt this in Figure]^ 
where we concentrate on one sightline, Z = 6 ± 1 degrees, 
for which both our simulation and APOGEE data display a 
peak at high wios- 

Panel E shows how the uios distribution varies when 
stars are sorted in four different age bins. Whereas the high 
uios population is very distinct for stars with ages below 1 
Gyr, this population becomes less distinct at higher ages. 
For stars with ages between 1 and 2 Gyr, it is still clearly 
visible, at 2-4 Gyr only a small feature remains and for older 
stars even this weak signature vanishes. The vanishing of 
the distinct high vios population goes hand in hand with 
broadening of the main uios peak indicating a heating of the 
population. 

To see how this evolves with simulation time, we plot 
in Panels A-D stars in these age bins at three different sim- 
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Figure 7. Line of sight velocity tiiog distributions N (vi^g) / Nmax of stars at longitudes Z = 14ibl, l = 6ibl and I = —6 it 1 degrees and 
latitudes |6| < 1 degrees. We assume a bar angle of 0 = 20 degrees, a solar radius of ij = 8.3kpc and use uiog bins which are all 8kms“^ 
wide. We compare four different types of histograms: blue considering all stars at i? < 4kpc, equally weighted (EW), red considering 
young stars with ages < 2 Gyr at ij < 4 kpc, equally weighted, green considering all stars along the sightline, equally weighted and black 
considering all stars along the sightline, weighted according to the age- and distance-dependent selection function SF(s,t) from Figure 
[^(SFW). The vertical dashed lines indicate the intervals for peak stars considered for Figure]^ 


ulation times, 7.3, 8.3 and 9.3 Gyr. We notice that there 
is much less evolution with simulation time than there is 
with age. The level of distinction of the high uios signature 
remains similar at all considered simulation times for each 
age bin. 

To strengthen these findings we plot in Panel F of Fig¬ 
ure]^ the population of stars born between 6.3 and 7.3 Gyr 
at three different times representing mean population ages 
0.5, 1.5 and 2.5 Gyr. This underlines that the evolution of 
the i>ios distribution with age shows the broadening of the 
main uios peak due to heating and the gradual loss of the 
high uioa peak. 

We again notice that the young populations of stars at 
9.3 Gyr show additional features at the opposite side of the 
main peak. Closer inspection reveals hints of this feature 
also at earlier times. Additionally, the main peaks move to 
higher uios. We will provide an explanation for these changes 
in Section 7. We note, however, that these additional peaks 
are only visible if we select young bar stars. When applying 
the selection function on all stars along the line of sight, 
as described in the next Section , these features, unlike the 
conventional high uios peaks, are too small to be detected. 


6 COMPARISON TO APOGEE DATA 


In the previous Section, we have shown that populations of 
young bar stars are expected to display high |uioa| features 


as observed by Nidever et al. (20121 and Babusiaux et al. 
( |2014[ ). In reality, not only stars with R < 4kpc will con¬ 
tribute to the APOGEE uios distributions, but also disc stars 
in front of or behind the bar. In the following, we directly 
compare the A-body simulation to the APOGEE data, ap¬ 
plying the selection function described in Section 2. 


6.1 The effect of the selection function 

We start with Figure where we compare the uios his¬ 
tograms for stars with galactocentric radii R < 4 kpc (see 
previous Section) to uios histograms obtained for the full 


sightlines. We do so for three longitudes, I = ±6 and 14 
degrees. The blue curves represent all stars within 4 kpc 
and the red curves show the subset of young stars, both are 
taken from Figure The velocity distributions for all stars 
along the sightline (equally weighted, green) are strikingly 
different. Their main uios peaks are significantly narrower 
and shifted by more than 50kms“'^ towards zero. This is 
due to the dominant contribution by fore- and background 
disc stars, which are kinematically cold and are observed 
with the local azimuthal direction nearly perpendicular to 
the line of sight. 

We noticed above that no distinct high |uios| feature is 
found for all stars at i? < 4 kpc. Interestingly this is not 
true for all stars along the line of sight. In this population, 
the feature is less distinct than the one for young bar stars, 
but lies at the same wios range. The histograms obtained 
when weighting all stars along the line of sight according to 
the selection function presented in Section 2 are shown as 
black lines. Overall, they resemble more the ones for all stars 
along the line of sight than to the ones for stars at i? < 4 kpc. 
When present, the high |uios| feature is strengthened by the 
selection function and the main peak is shifted to slightly 
larger |uios|. 

From these comparisons we conclude that the Uios distri¬ 
bution is not only shaped by stars in the bar region, but by 
all stars along the line of sight. This is illustrated by the se¬ 
lection function weighted (SFW) distribution of stars along 
the line-of-sight c{s) as a function of distance from the sun 
s in the left panels of Figure (dashed lines). Clearly, the 
SFW distribution functions do not simply reflect the den¬ 
sity distributions along the pencil beams, but are shaped 
by the age and distance dependent selection function modu¬ 
lated by the beam area increasing as s^. Of the three shown 
sightlines, only the one which crosses the bar most centrally 
(Z = 6 degrees) has its absolute maximum coincident with 
the bar distance. The loci of young stars, the highest densi¬ 
ties of which are by construction at a galactocentric radius 
R ~ 4 kpc and thus at distances s ~ 4 and 12 kpc are visible 
as local maxima. Moreover, spiral structure creates minor 
peaks. 
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Figure 8. Left panels: The selection function weighted distribution of all stars in a sightline (dashed lines) and of the stars contributing 
to the high |uios| peaks (solid lines) as a function of their distance from the Sun s, c(s)/cniax- We show longitudes Z = 14 it 1 (red), 
Z = 6 it 1 (orange) and Z = —6 it 1 (blue) degrees. High |uios| windows are defined as 180 < uioa/kms”^ < 250 for Z = 14 it 1, 
220 < uioa/kms“^ < 280 for Z = 6 it 1 and —300 < uioa/kms“^ < —200 for Z = —6 it 1, as indicated in Figure [3 Right panel: The 
location of the high |uios| peak stars in the x — y plane. We show the general geometry of single sightlines from the Sun (yellow star) into 
the bar (stellar mass density contours from Figurej^in black). Arrows indicate galactic and bar rotation. The black circle marks R = 4 kpc. 
The magenta areas show the distance areas which contain 95 per cent of the SFW stars contributing to the high |uios| peak windows for 
stars moving away from the sun, cyan areas are for high |uios| stars moving towards the sun at negative longitudes. Darker shaded areas 
show the areas of ±0.5kpc centred on the peaks of c(s). Note the position of the fast moving stars with 180 < uioa/kms“^ < 250 at 
Z = 14, which indicates the increasing contribution of stars outside the bar. 


The high nios feature appears connected to young bar 
stars. To verify this assumption, we plot in the left panels 
of Figure the contribution to the high |uios| features of 
SFW stars for the three sightlines of Figure [T] as a function 
of distance from the Sun s (solid lines). The intervals we 
use as definitions for these features are indicated by dotted 
vertical lines in Figurej^ For I = ±6 degrees, we were guided 
by the distinctly visible feature, for I = 14 degrees, where 
we find no feature, we use the estimate by |Nidever et al.| 

| M2| ). 

At each sightline in Figure the regions of contribu¬ 
tions to the high |uios| stars are much narrower than the 
SFW distribution of all stars in the pencil beam. Moreover, 
the maxima of contribution to the high |uios| stars are (apart 
from I = 6 degrees) offset from the maxima of the SFW dis¬ 
tribution of all stars. The distance distribution of contribut¬ 
ing stars c(s) is narrowest for I — 6 degrees and widest for 
I = 14 degrees. The main peaks are at s = 7 — 8.5 kpc on 
the rear side of the bar (front side for I = —6 degrees). To 
visualise where these contributions arise in relation to the 
bar, we use the right panel of Figure]^ where we show the 
position of the Sun in yellow, the considered sightlines in 
blue, red and orange and the surface density contours of the 
stellar mass in the bar region at 7.3 Gyr in black (see Fig¬ 
ure]^. The 95 per cent contribution regions of high |uios| 
stars for each sightline are shown in magenta for stars mov¬ 
ing away from the sun and in cyan for stars moving towards 


the sun at negative longitudes. We highlight the regions of 
±0.5 kpc around the peak distances in darker shades. For 
I = ±6 degrees, these regions are fully within the bar region 
and strictly on one side of the bar (rear side for I > 0, front 
for Z < 0), as the fast stars travelling alongside the bar have 
to lose angular momentum on the way in and gain it on the 
way out. The contribution region is narrower for I = 6 de¬ 
grees as this sightline crosses the bar at a greater angle than 
the one for Z = — 6 degrees. 


This indicates that the high |uios| features found by 


Nidever et al. (20121 correspond to stars which are asso¬ 


ciated with the bar region and not with disc orbits. This 


confirms the ideas of Nidever et al. (20121 and Molloy et 


al. (20151 and rejects the notion of Li et al. (20141 that the 


observed features are not related to the bar. 


At Z = 14 degrees there is a tentative feature in observed 
data, but not in our simulations. If we select the highest uios 
stars in the simulation we find that the main contributions 
come from the edge of the bar, but extend far into the disc 
region behind the bar. At these longitudes, the inner parts of 
the disc align with the line of sight and start to contribute to 
the high uios window. Should the high uios feature at Z = 14 
degrees prove to be real, a clear attribution may only be 
possible with detailed analytical models. 
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Figure 9. A comparison with APOGEE data. For 12 sightlines for which data from APOGEE are available we compare line of sight 
velocity 'Uiog distributions N{viag) to the data discussed in Section two. The APOGEE data shows the actual number of stars per bin, 
the simulation data are normalised to yield the same total number. The APOGEE -uiog bins are 25kms“^ wide, whereas the simulation 
bins are 8kms“^ wide. For comparison to APOGEE data (magenta), we show simulation histograms for equally weighted (EW, green) 
and selection function weighted (SEW, black) stars along the sightline. We assume a bar angle of </> = 20 degrees and a solar radius of 
R = 8.3 kpc. 
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Figure 10. Line of sight velocity vi^s distributions Ai(uios)/7Vmax as seen from R = 8.3 kpc at a bar angle of </< = 20 degrees at a 
simulation time of 7.3 Gyr. ujog bins used here are all 8kms“^ wide. We consider sightlines centred on I = ±6 degrees, which are both 
2 degrees wide in 1. In each panel we consider four latitude bins: |b| < 1 degrees (black), 1-2 degrees (blue), 2—4 degrees and > 4 degrees 
(red). Left panels: Here we consider histograms of all stars along the sightline, weighted by the selection function depicted in Figure]^ 
Right panels: Here we consider equally weighted, young stars (ages < 2 Gyr) at i? < 4 kpc. 


6.2 Detailed comparison with APOGEE 

We now directly compare our model with APOGEE data. 
Remember that we have done nothing at all to fit this spe¬ 
cific model galaxy to the MW. There may be significant 
differences in length, mass, age and detailed history of the 
bar etc. Hence the goal is to qualitatively identify features 
and not to do a full quantitative analysis. 

In Figure we compare vios histograms for all stars 
along the sightlines, equally weighted (EW, green) or se¬ 
lection function weighted (SFW, black) to the APOGEE 
data described in Section 2 for various sightlines. We dis¬ 
play in-plane data at |&| < 1 degrees in the upper two rows, 
and sightlines at higher |6| in the bottom two rows. Despite 
its limitations, the simulation fits the data surprisingly well. 
The SFW models provide a better fit to APOGEE data than 
the EW models. Their high uios features are more distinct 
and the EW histograms are offset from the data at large I 
as well as at large b. 

The sightlines with I > 14 degrees are increasingly dom¬ 
inated by disc stars. There, the parts aligned with the line of 
sight would bury any distinct velocity feature from the bar. 
In concordance with this expectation, there is no discernible 
feature in the corresponding in-plane APOGEE data, and 
also not in the model. In contrast, at Z = 3 — 10 degrees, the 
high velocity feature is prominent both in the MW data and 
our simulation. Taking into account the significant Poisson 
noise in the observations, there is little difference at Z ^ 6 
degrees. At Z = 10 degrees, the feature is less distinct in the 
simulation. Interestingly, at Z = 16 degrees the SFW his¬ 
togram provides a worse fit than the EW histogram, as it 
misses the position of the main peak, which at other longi¬ 
tudes is well reproduced. 

When we look slightly above and below the plane at 
l&l = 1 — 3 degrees, asymmetries at high uios are present in 
APOGEE data at Z = 5—15 degrees, but again not at greater 
1. The tentative Z = 14 peak in the data was highlighted by 
), but has a very weak counterpart in 
SFW model histogram here drops less 
steeply than the data on the high velocity side of the main 
peak, as is the case for Z = 16 degrees (mentioned above). 


Nidever et al. (2012 


the simulation. The 


These are the strongest discrepancies between model and 
observations. There may be various explanations for this: 

i) We know that the Galactic potential will not be per¬ 
fectly matched. In particular for disc contributions, the loca¬ 
tion of the main uios peak is determined by the radial poten¬ 
tial gradient (in idealisation the circular speed). If this dif¬ 
fers, simulation and data will show a systematic difference. 

ii) We know that the structure in particular in this region 
will show differences between simulation and observed data, 
e.g. caused by a difference in bar length. A longer bar would 
mean that we see regions further away from the bar tips. We 
can achieve such an effect also by lowering the solar radius. 
Indeed, if we do so, the asymmetry becomes more prominent. 

iii) Another explanation could be the dense ring of molecu¬ 
lar clouds at i? = 4kpc in the MW ( Dame et al.||1987 |. This 
ring is associated with an intense radial peak in star forma¬ 
tion, which is not matched by our simulations. The selection 
function of APOGEE is very sensitive to young stars. This 
feature will have no strong consequence when the sightline is 
perpendicular to circular orbits at R ~ 4 kpc (at small Z) or 
outside the region (Z > 30), but sightlines at Z ~ 16 degrees 
will be strongly affected. As they pass through R — 4 kpc 
at an intermediate angle to the local azimuthal direction, 
the enhancement of these stars will drive up stellar densi¬ 
ties at intermediate Uios and shift the main peak. Overall, 
the complicated contribution function of stars at these sight¬ 
lines prevents any strong conclusions to be drawn from this 
discrepancy. 

Moving even further away from the plane, we consider 
latitudes |6| is 3 degrees in the lower row of Figure]^ There 
is no discernible feature, either in observations or simulation. 
The slight asymmetry is linked to the rotation of the bar. 
So our model provides reasonable fits to the data in and out 
of the plane. 


6.3 Varying the parameters 

We now examine the behaviour at different latitudes. Figure 
|10| shows simulation data for two sightlines Z = ±6 degrees 
and varying latitude. We do this both for SFW histograms 
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Figure 11. Normalised line of sight velocity tiiog distributions n(iiiog) of stars at 7.3 Gyr simulation time and latitudes |f)| < 1 degrees 
assuming a solar radius iJsun = 8.3 kpc. If available APOGEE data are overplotted in magenta. In each panel we consider bar angles (p 
between 10 and 55 degrees. Upper row. All stars in the sightline, weighted by the selection function presented in Figure^ Lower row. All 
stars within i? < 4 kpc and with ages < 2 Gyr, equally weighted. Left column: Longitude 1 = 10 it 1 degrees. Middle column: Longitude 
1 = 6 it 1 degrees. Right column: Longitude I = —6 it 1 degrees. 


for all stars along the lines of sight (left panels) and for 
EW young central stars (right panels). The difference be¬ 
tween these two types of histograms is even more striking at 
higher |6| where for young central stars they are significantly 
broader. This indicates that young bar stars play a small role 
in observations at higher latitude. Note also that the number 
of young bar stars decreases strongly with increasing lati¬ 
tude in the simulations, which explains the noisier nature of 
the corresponding histograms. We find that the exponential 
scaleheight for stars with ages < 2 Gyr, which contribute 
to the high uios features is hz ~ 70 — 80 pc, compared to 
hz ~ 180 pc for all stars at the ends of the bar. 


All panels show that distinct high |wios| features are only 
present close to the plane. At |6| = 1 — 2 degrees a mild 
feature is still present, while at greater b there is nothing left. 
This explains why Zoccali et al. (20141 did not discover such 
features in their data, which are mostly at |&| ^ 2 degrees. 
Their only field at smaller latitude is at 1 = 0, where we do 
not expect a peak anyway. This can be explained by the fact 
that the high uios populations are dynamically cool and thus 
also live at small vertical distances from the galactic plane. 


So far, we have assumed a bar angle of ((> = 20 degrees. 
In the literature, estimates range from 15 to 40 degrees. |Mol-| 


loy et al. (20151 recently claimed that the high uios features 


observed by Nidever et al. (20121 could help to constrain (p 
and would favour low values 0 ~ 15 degrees. In Figure m 
we thus attempt to study the influence of (p on uios distri¬ 


butions. We consider angles between 10 and 55 degrees for 
three sightlines and both SFW stars (upper row) and young 
stars with R < 4kpc (lower row). 

We see that both for positive and negative Z, the high 
|uios| features are less distinct for large (p. This can be ex¬ 
plained by the larger angle between the sightline and the bar. 
As the high velocity stars move along the bar, their projected 
Wios decreases. At larger (p, stars on different orbits moving 
at greater angles to the bar can be seen at higher uios. As 
they are less numerous, the feature becomes less prominent. 
This is especially clear for young stars at i? < 4 kpc. If one 
tried to fit multiple Gaussians to the Uios distribution, as did 


Nidever et al. (20121 and Molloy et al. (20151, one would 


find that the high uios peak becomes broader and centred at 
lower Vies as reported by Molloy et al. (20151. 

In the upper panels of Figure [m APOGEE data are 
plotted over the SFW model. We should keep in mind the 
uncertainties in the observations and the fact that the struc¬ 
ture and age distribution of the model bar differ in an un¬ 
known extent from the MW bar. We thus conclude that the 
APOGEE data put only mild constraints on the bar angle. 
They favour small angles, but only strongly disfavour angles 
(() > 40 degrees. 

The estimate by Molloy et al. (20151 is also intrinsically 
overly confident, as they ascribe the high wios feature to one 
orbit family only, whereas, as we will show in the next Sec¬ 
tion, several families can contribute. Moreover, degeneracies 
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Figure 12. Left panel. Line of sight velocity tiiog distributions N(vi^s)/Nmax of young stars at ij < 4kpc as seen from R = 8.3 kpc 
with a bar angle 0 = 20 degrees at longitudes Z = 6 it 1 degrees and latitudes |6| < 1 degrees. The stars were selected as members of 
the high velocity peaks, defined by 210 < uios/kms“^ < 280 at 7.30-7.35 Gyr (blue) and by 230 < uios/kms“^ < 300 at 9.30-9.35 Gyr 
(green). Then they were followed for 500 Myr with a snapshot separation of 1 Myr and each time they were in the sightline Z = 6 it 1, 
their uiog was recorded to produce the shown uiog distributions. Right panel. The positions of stars at 9.30-9.35 Gyr seen at longitudes 
Z = 6 it 1 degrees and latitudes |Z)| < 1 degrees and either in the high uiog peak (230 < uios/kms“^ < 300, black dots) or the low uiog 
peak (—150 < uios/ kms“^ < —50, red dots). Each star is recorded only once, on its first detection in the corresponding velocity window. 
The sun is at x = —8.3 kpc and y = 0 here. Overplotted in green is an orbit over 500 Myr, which during this time alternately contributes 
to both peaks. The arrows indicate the sense of motion. 


with the potential, the pattern speed, the bar structure etc 
complicate their determination of 0. 


7 TRACKING ORBITS 

To understand better the origin of the high and low uios 
peaks discussed above, we focus in this section on the orbits 
producing these features. 


7.1 High and low uios features 

We created snapshots of the simulations at a frequency of 1 
per 1 Myr from simulation time 7.3 Gyr onwards. We first 
consider the following exercise: What would the uios distri¬ 
bution look like if it were only made up of stars that at 
some point contribute to the high uios peak. We therefore 
select all stars with ages < 2 Gyr and R < 4 kpc that at 
7.30-7.35 Gyr contribute to the high uioa peak at Z = 6 ± 1 
degrees and |fo| < 1 degrees. We then follow these particles 
for a further 500 Myr and construct a uios distribution at 
the same I and b ranges only from these stars. We do the 
same 2 Gyr later for a new set of stars selected in the same 
way. Figures and have already shown that the maxi¬ 
mum (minimum) uios observed increase (decrease) by about 
20kms“^ within these 2 Gyr. This is because of a growth of 
both the potential well and the bar and also because of the 
slowdown of the bar which allows orbits to have lower angu¬ 
lar momentum and reach positions closer to the centre and 
thus cover a larger potential difference. Gonsequently, we in¬ 
crease the high uios window from 210 < Uioa/kms“^ < 280 
to 230 < uios/kms“^ < 300. 


These uioa distributions are shown in Figure Unsur¬ 
prisingly they have their main peaks at the selection po¬ 
sitions. Both distributions clearly show a distinct peak at 
negative uios made up by orbits returning along the other 
side of the bar. This peak was already visible in Figure 
for 9.3 Gyr. The increase in maximum absolute velocities 
between 7.3 and 9.3 Gyr widens the gap between the peaks 
and shifts the main peak of the overall distributions for stars 
within R < 4 kpc to higher uiog. Gonsequently the negative 
uios peak at 9.3 Gyr for Figure [^becomes visible. The neg¬ 
ative peak comprises approaching stars while the positive 
peak features stars moving away from the sun. Figure 
reveals a third peak for 7.3 Gyr indicating that a group 
of orbits has a more complicated structure than alternately 
contributing to the two peaks. 

In the right panel of Figure we visualise the origin 
of the high and low vios peaks at 230 < Uios/kms“^ < 300 
and —150 < uios/kms“^ < —50 for young bar stars at 9.3 
Gyr from spatially separate components along the line of 
sight. Here we plot the spatial positions of all stars in the 
positive uios peak at 9.3 Gyr in black and of all stars in the 
negative uios peak at 9.3 Gyr in red. The components clearly 
separate to the expected sides. One orbit family that would 
contribute to both peaks are the xi orbits. In the rotating 
frame of the bar, these orbits move along opposing sides of 
the bar major axis and turn around at the bar tips. They 
would obviously produce the two peak feature in the left 
panel of Figure |12[ as was e.g. discussed by [Molloy et al.| 
( |2015[ ). We overplot in green one such orbit extracted from 
the simulation which alternately contributes to both peaks. 
However, as we will show below, the high uios features in our 
simulation are made up of a more complicated population 
of orbits from different orbit families. 
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Figure 13. For this Figure, we track stars which at 7.30-7.35 Gyr are younger than 2 Gyr and at galactocentric radii R < 4kpc, and 
contribute to the high peak (210 < uioa/kms“^ < 280) at / = 6 it 1 and |fe| < 1 degrees. We track them for 2.4 Gyr afterwards. We 
record the first, second and third time (ri, T 2 and T 3 ) at which they are at t = 6 it 1 and 210 < uioa/kms“^ < 280 and use it to calculate 
the first and second return time intervals t\ = T 2 — ri and t 2 = — T 2 . Panel A: The histogram for ti. The peak selection windows are 

marked by vertical lines in red, green, blue, purple and orange. Panels B-F: t 2 histograms for orbits selected to be in one of the ti peak 
windows marked in Panel A. The colour of each panel corresponds to the colour of the selection window. The selection windows from 
Panel A are marked by dotted vertical lines. 


7.2 A simple classification for high uios orbits 


To get a better understanding for the orbits that contribute 
to the peak, we track all orbits that contribute to the high 
Wios peak at Z = 6 ± 1 and |fe| < 1 degrees at 7.30-7.35 Gyr 
defined by 210 < uios/kms”^ < 280 for young stars with 
ages < 2 Gyr and radii R < 4kpc. If an orbit is similar to a 
closed orbit, we will expect it to return to the same region in 
phase space periodically. Our selection of sightline plus uioa 
range can be seen as a rough proxy for a phase space region. 
Considering that the bar grows and the potential changes, 
a more sophisticated definition is complicated. Therefore we 
record the time interval needed by an orbit to return to 
the sightline and uios window. The histogram for the first 
return time ti is plotted in the Panel A of Figure |13| We 
see that this distribution shows three distinct peaks at circa 
110 (red), 190 (green) and 280 Myr (blue), which contain 
approx. 23, 13 and 28 per cent of the stars. There are also 
two broader peaks at 450-650 Myr (14 per cent, purple) and 
650-800 Myr (7 per cent, orange). Note that ~ 13 per cent 
of orbits never return or have return times > 800 Myr. Were 


all these orbits simple x\ orbits, we would expect a simpler 
structure. 

In Panels B-F of Figure we look at histograms of the 
second return time t 2 for stars selected to have first return 
times ti in the five peaks found in the ti histogram. To some 
degree, all these t 2 histogram resemble the t\ histogram, as 
they reproduce one to three of the previously found three 
narrow peaks. Above 450 Myr, the t 2 histograms all seem 
dominated by noise, but all of them contain a similar fraction 
of such orbits, varying between 25 and 31 per cent. A first 
clear conclusion is, that for each histogram, the fraction of 
orbits that remain in the same return time peak is below 
50 per cent. This shows that most orbits do not show a 
steady periodicity in our definition, but rather combine a 
low number of possible return times. 

To connect the shape of these histograms to types of 
orbits, we extracted random sets of orbits with similar com¬ 
binations of return times ti and ± 2 , for example contributors 
to the green peak for ti and to the blue peak for t 2 - In ad¬ 
dition we extracted a set of orbits that do not contribute to 
return time peaks. From these we select by eye 42 typical 
orbits. They are shown in Figurej^ Note that due to differ- 
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Figure 14. 42 Orbits in the x — y plane of particles that at 7.3 Gyr contribute to the high uios peak at / = 6 di 1 degrees. The orbits 
are shown from when they first leave the sightline until the second time they return to the sightline and contribute to the high peak 
again. The numbers in the corners of the panel indicate the times between contributions to the peak in Myr. If these fall into one of the 
return time peaks identified in Figurethey are coloured accordingly, if they don’t, they are coloured black. The orbit colour coding 
indicates time, changing from black via light and dark blue to red with increasing time. 
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ent prominence of the peaks this figure does not represent 
orbit types at their underlying contribution fractions. 

The stars with shortest return times ti (red peak) show 
a high fraction (45 per cent) of stars with ti ~ t 2 . Addition¬ 
ally, 30 per cent of the stars have t 2 in the area of the green 
or blue peak, but these peaks are not separately visible. The 
first row of Figure reveals that these stars are essentially 
all on like orbits. If they have t 2 ^ 1 1 , it is mostly because 
the orbits are not exactly closed. They can miss the high uios 
window as their directions of motion and velocity vary when 
they cross the sightline. Moreover, A^-body noise can cause 
the integrals of motion of the particles to fluctuate. 

Stars from the second (green) ti peak show three t 2 
peaks at the locations of the three narrow ti peaks which 
are equally populated (~ 25 per cent each). Most stars are 
again on xi-like orbits, but there are also the first objects 
from more complex closed orbit families with Pretzel-like 
structures j Portail et al.|2015b l, such as the 202-213 or 216- 
229 objects. At higher return times, Xi orbits are a minority. 
If they have such return long return times, it is because they 
miss uios due to the fluctuations mentioned above more often 
than not, e.g. because they have lower energies and only in 
exceptional cases reach high enough uios. 

Pretzel-like morphologies are the most frequent pattern 
found at longer return times. For blue ti stars, there is again 
a high fraction (46 per cent) of ti ~ t 2 stars, which tend to 
shows such patterns, as seen in the fourth row of Figure 
|14[ Generally, the structures become more complicated and 
more irregular as ti increases. Clearly some orbits show 3:1 
characteristics with a changing orientation (objects 421-402 
or 506-862), but there are also more complicated, but appar¬ 
ently regular structures, such as the object 273-176, which 
was labelled ’Q’ in Voglis et al. (20071, and objects like 176- 
324 or 424-103 which also combine 3:1 and x\ features. Other 
orbits appear to change structure with time, an indication 
of chaotic behaviour. 

If we combine the results from above, we estimate that 
the fraction of stars in the high uios peak that are on xi 
orbits at the given time is between 40 and 50 per cent. If 
we determined the full population of stars, which over a 
longer timescale (~ 1 Gyr) periodically contributed to the 
high uios peak, the fraction of X\ stars would be even lower, 
as they have the shortest return times. This invalidates the 
original assumption of Molloy et al. (20151 that the high uios 
feature is only caused by 2:1 orbits. There are other major 
contributors, most prominently the Pretzel-like orbits, but 
also 3:1 or more complex orbits. We note that |Molloy et ah] 
| |2015[ ) have recently updated their paper in response to our 
work and now also find that higher order orbit families (they 
name 3:1 and 5:2) can make significant contributions to the 
high uios features. 

We note that we have attempted to analyse the vertical 
structure of orbits in connection to the in-plane one and ex¬ 
perimented with frequency analysis, but concluded that the 
results were to diverse or noisy to help with the understand¬ 
ing of the high uios features. 


T.3 Are orbits evolving away from high uios orbits? 

In our discussion of Figure |13[ we already noted that around 
10 per cent of stars which contribute to the high uios peak at 
a given time and sightline, need 1 Gyr or longer to achieve 



Figure 15. For this Figure, we track all orbits from the his¬ 
togram in Panel A of Figure |13| We exclude orbits that never 
return to the high ujoa window or have very long return times 
ti > 0.8 Gyr. We record each passage of the sightline I = 
6 ± 1 degrees (total number A^indow(t)) and each contribu¬ 
tion to the high uioa peak 210 < t)ios/kms“^ < 280) during 
such a passage (total number Afpeak(t))- We plot the fraction 

^(t) — Apeak (t)/Awindow(t). 


this again. The t 2 histograms have also shown that ~25 per 
cent of stars that initially had return times < 400 Myr need 
more than 500 Myr to return again. So orbit populations 
slowly change in a way that reduces the proportion of their 
time during which they contribute to the high uioa peak. 

To understand why the high uios peak vanishes with age, 
we conducted the following exercise: We tracked all orbits 
from the histogram in Panel A of Figure As the orbits 
that never return or have very long return times (defined 
as > 0.8 Gyr) are likely random contributors, we exclude 
them. We record each passage of the sightline (total number 
Awindow(f)) and if during this passage the orbit contributes 
to the high wioa peak (total number Apeak (f))- The fraction 
f = Apeak (t)/Awindow(t) gives a measure of how much the 
population of orbits is contributing to the peak at time t. 
A priori one would expect a ratio v ~ 0.5 for xi orbits as 
they are approximately axisymmetric along the bar passing 
once through the sightline while approaching and once while 
moving away and contributing to the high vioa peak. Con¬ 
sequently, smaller values for v are expected for more com¬ 
plex orbits. This is slightly complicated by the geometry and 
the slower absolute velocity of stars while approaching: The 
front side of the bar is geometrically disfavoured, as the pen¬ 
cil beam covers a larger area at the rear side. On the other 
hand, the approaching stars are closer to the tip and hence 
slower, taking more weight in the observations. Moreover, it 
is also possible that stars moving along the rear side have 
high enough uios only for a fraction of the time they spend 
in the sightline. If we select xi orbits which at every second 
crossing of the sightline contribute to the high uios peak, we 
empirically find a mean value of iz « 0.45, close to the simple 
estimate. 

From the return time analysis shown above, we would 
expect 1 / < 0.45, as most orbits exhibit a more complex 
behaviour. iz(t) is shown in Figure [T^ Initially, i/ fluctuates 
periodically. This can be understood as most orbits return 
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Figure 16. Orbits of captured particles in the x — y plane. These stellar particles are born between 7.3 and 8.3 Gyr and contribute 
to the high uiog > 220kms“^ feature at Z = 6 it 1 degrees and 8.3 Gyr. The numbers in the top left corners indicate the time in Myr 
between birth and first detection in the high uiog window at Z = 6 it 1 degrees. The colour coding indicates time, from black at birth via 
light and dark blue to red at first detection. 


to the high vios peak after specific characteristic times (see 
Figure [l3| ) and they have been selected at correlated phases. 
The average value of u over the covered 2 Gyr drops from 
ly ~ 0.22 to V ~ 0.12. This shows that the fraction of orbits 
contributing to the high vios peak continuously decreases. 

Most of this decrease is caused by diffusion of stars away 
from the contributing orbits, but it is not clear, if our simula¬ 
tion handles this diffusion correctly. Sources for this diffusion 
are chaotic orbit behaviour, two-body noise and fluctuations 
in bar parameters. Our resolution is likely not sufficient to 
cover chaotic orbital behaviour correctly. For the two-body 
noise, it is unknown if the level of Ai-body noise is realistic 
compared to the MW (bearing e.g. Giant Molecular Glouds 
near the bar tips, and dark matter substructure). In addi¬ 
tion, our simulation will not have sources for bar fluctuations 
like tidal torques exerted by satellites. 

The important point to take away is that the vanishing 
of the high uios peak in older stars is not only due to the 
broadening of the main peak, but also due to active loss of 
stars from high uioa orbits. 

T.4 Orbit capture by the bar 

By design, new-born star particles in our simulation are re¬ 
leased on near circular orbits, where the circular speed was 
determined to balance the local radial component of the 
gravitational force. Obviously, the young stars contributing 
to the high uioa peak discussed above are not on circular 
orbits. Just as obviously, there are no stable circular orbits 
near the bar. This is why we set an inner radius limit for new 
born stars, so that no new star particles are inserted within 
the main bar region. This cutoff is, of course, to some de¬ 


gree arbitrary and young stars born close to the bar can be 
captured by the bar, resulting in the young bar populations 
studied in the previous Sections. Bars in real galaxies can 
also grow by capturing stars and thus changing their orbits. 
So the young stars that end up in the high uios peak are 
predominantly the result of continuous bar growth. 

In Figure |16[ we show six stars that at 8.3 Gyr con¬ 
tribute to the high uioa > 220kms“^ feature and were born 
within 1 Gyr before that time. We find that the shortest 
time it takes to appear in the high wioa peak is similar to 
the shortest return time peak, 110 Myr. The capture times 
are rather equally distributed between this lower limit and 
the maximum 1 Gyr due to the selection criteria. The shown 
star particles were selected at different capture time inter¬ 
vals. They start out on non-bar orbits and at some point 
transition to xi-like or Pretzel-like orbits. 

Clearly stars that end up in the high uios features are 
born at radii that correspond to the tip of the bar, so they 
belong to the innermost inserted particles. To visualise this 
point, we plot in Fignre the distribution of birth radii 
I?birth of high uios stars. For this Figure we select stars which 
were born after 7.3 Gyr and contribute to the high wios fea¬ 
ture at Z = 6 degrees at 8.3-8.5 Gyr and/or at 9.3-9.5 Gyr. 
The union of the samples found at these two time intervals 
comprises ~ 11000 stars, for which we determine birth radii. 
The histogram for all stars is shown in blue, the one for the 
subsample of stars born after 8.5 Gyr is shown in red. We 
notice that the majority of stars is born just outside the bar 
at i? ~ 4kpc. As seen in Figure the bar radius fluctu¬ 
ates mildly, which allows smaller birth radii. The subsample 
of stars which were born at later times has its lower cutoff 
slightly further out due to the continuing growth of the bar. 
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scale height remain small. Consistently, both the MW and 
our simulation display a discernible high vios feature only at 
observed latitudes b <2 degrees. 

Above, we have learned that for populations of stars 
which at a given time contribute to the high uios feature, 
the contribution fraction slowly declines with time. As the 
simulation evolves, the high uioa feature, however, remains 
detectable. This is because of ongoing star formation close 
to the bar region and the continuous capture of young stars 
on Xi-like or Pretzel-like orbits. The young populations on 
such orbits then slowly evolve towards more complex orbit 
populations with a decreasing fraction of high uios stars. 

We note that another process which could produce 
young stars on high uios orbits is not included in our model: 
In real galaxies, molecular clouds can exist on bar orbits 
i Sheth et al.||20d2 |. Stars born in these clouds would make 
up kinematically cool populations born directly on high ve¬ 
locity bar orbits and could thus contribute to the high uios 
features. 


Figure 17. For this Figure we select stars which were born after 
7.3 Gyr and at 8.3-8.5 or at 9.3-9.5 Gyr contribute to the high vi^s 
window at I = 6 degrees. The union of the samples found at these 
two time intervals comprises ~ 11000 stars. The distribution of 
birth radii 7?birth of all ~ 11000 stars is shown in blue. The 
histogram of stars born after 8.5 Gyr and contributing at 9.3-9.5 
Gyr is shown in red. The vertical dashed lines mark the positions 
of the inner 4:1 (ultra-harmonic) resonance at 7.3 and 9.3 Gyr. 


The outer tail is wider for older stars, indicating that stars 
which are older at capture can have larger birth radii. This 
is because they had more time to migrate radially inwards. 


Resonances can play a role in capturing orbits (e.g. Bin 
ney &: Tremaine|2008 1. The corotation resonance of our bar 


at the analysed time intervals is at J? = 7 — 9 kpc, but no 
stars from these regions are captured. We find that the ultra¬ 
harmonic resonance (UHR), 4(f2 — flbar) = in our model 
lies just outside the bar. This is likely not coincidental, as the 
UHR can limit the radial extent of stable xi orbits and thus 
the size of the bar ( Patsis et al.||1997 |. The resonance may 
play a role in shaping the spatial distribution from which 
stars are captured by the bar, as while the bar grows, the 
resonance moves outwards, and stars on orbits just outside 
the bar get scattered, absorbed or released by it. These pro¬ 
cesses can transport orbits to the regions of phase space oc¬ 
cupied by high velocity bar orbits. However, the exact role 
of the resonance in the capture process is unknown, as the 
bar is strong, evolving and the potential is not smooth. The 
capture of orbits can be understood in simpler terms (simi¬ 
lar to Lynden-Bell 19791, in that stars come in at the right 
angle to the bar to lose angular momentum and get trapped. 
A detailed study of capture processes is beyond the scope of 
this work. 

Capture of stars onto these orbits along the bar is 
mostly drawn from kinematically cool populations that are 
close to the Galactic plane, which have a small random en¬ 
ergy compared to the potential disturbance and are at low 
altitude, where the disturbance is strongest. While these 
stars drastically change the in-plane geometry of their orbits 
during capture by the bar, their vertical energy and hence 


8 CONCLUSIONS 


In this paper, we have studied how an N-body-l-SPH model 
of a growing Milky Way (MW) like barred disc galaxy com¬ 
pares to the MW in terms of line of sight velocities uios of 
stars seen from the earth and in the direction of the bar. The 
model uses an approach that has not yet been used in this 
context. We start with a light (one tenth of the final mass) 
and short scale-length disc of gas and stars in equilibrium 
within a live dark halo. This disc is grown by continuously 
adding star and gas particles on near circular orbits in a con¬ 
trolled way. The mass and size growth histories are matched 
to expectations for the MW, and indeed at simulation time 
7-10 Gyr the simulated bar is similar in structure to the MW 
bar. 


Our simple prescription for the evolution of the galaxy 
in our N-body simulation lowers the computational cost and 
also circumvents the uncertainties in the parametrisation 
of sub-grid model approaches to star formation and stellar 
feedback. The general dynamics of gas and stars are ex¬ 
pected to be reasonably represented, but we do not cover 
varying gas fractions, supernova induced turbulence and ef¬ 
fects of multiple gas phases. The simulation does not pro¬ 
duce any sensible interplay between gas and stars: there is 
e.g. no enhanced star formation at places with high gas den¬ 
sities. External galaxies show indications for star formation 
associated with the gas flows along the edges of a bar (the 
equivalent of the MW’s 3 kpc arms; Sheth et al. 2002 I, which 
depend strongly on the dynamical state of the bar (see e.g. 
Friedli fc Benz|1995l. However, these processes are not cov¬ 


ered in our simulation. On the positive side, our approach 
allows us to set a reasonable star formation history, to keep 
spurious sources of heating under control, and hence to ob¬ 
serve the dynamical effects of bar formation with less scatter. 

We use this simulation to interpret line of sight veloc¬ 
ity (uios) data for the central and inner disc regions of the 
MW from the APOGEE survey. To enable a robust com¬ 
parison, we apply a simple population synthesis to model 
the APOGEE selection function. Apart from minor discrep¬ 
ancies at a minority of sightlines, the model fits the data 
surprisingly well without any scaling or adjustment. The 
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simulation not only reproduces the positions and dispersions 
of the main peaks in the uios distributions at each sightline, 
but also displays the high uios features found by|Nidever et| 


al. (20121 within the central bar region. 


Our simulation clearly identifies these high uios features 
with stars travelling near the plane alongside the bar. By 
showing that these are kinematically cool features with pref¬ 
erentially young stars in the plane, we naturally explain why 
the high uios structure was seen by Nidever et al. (2012 \ and 
|Babusiaux et al.|(|2014|) in observation fields at |6| < 2, but is 
absent in the higher-latitude fields studied by |Zoccali et al.| 
(20141. By the strong preference of the APOGEE selection 
function for young stars, the high uios features get addi¬ 
tionally enhanced and have comparable strength and shape 
in our simulation and the APOGEE observations. We also 
determine at which velocities low uioa peaks at negative I, 
produced by stars moving towards the Sun along the front 
side of the bar, should be visible. So far, only hints of these 
features have been seen by Babusiaux et al. (20141. 

We show that, unlike for the high |uios| features, the po¬ 
sition and width of the main peaks of the uios distributions 
are dominated by fore- and background stars, which are not 
associated with the bar. We also study a variation of our pa¬ 
rameters, in particular the bar angle (p. Other than |MolIoy| 
|et al.| ( |2015[ ) we feel unable to put strong constraints on the 
observed bar angle: while large ij) ^ 40 deg are clearly dis¬ 
favoured as the majority of high velocity stars get a too large 
projection angle with the line of sight, no reliable constraint 
can be found for smaller values of rp, since change in veloci¬ 
ties is degenerate with our general uncertainty about the bar 
length, pattern speed and strength. Additionally, we show 
that the high uios feature is detectable at simulation times 
between 7 and 10 Gyr and is thus not connected to specific 
events. The feature changes mildly over this period due to a 
growth of the bar in length and mass, but the change is too 
small to allow any conclusions about the properties of the 
MW bar. 

We note that on two sightlines, namely I = 14 and I = 
16 degrees, the simulation shows a mild difference to the 
observations, missing the tentative high uioa peak at i = 14 
and predicting the main peak of the velocity distribution 
at too high velocities at both sightlines. Possible origins for 
the discrepancies are a longer bar length for the MW or the 
4kpc molecular ring in the MW. A decisive answer on this 
question would require detailed analytical models that can 
cover the full space of possible parameters. 

A detailed analysis of the stars participating in the high 
wios features in our model reveals that (contrary to the origi¬ 
nal view of |Molloy et al.|2015[ ) only about half of the contri¬ 
bution is made by stars on xi orbits, with additional large 
contributions from Pretzel type (Portail et al. 2015b I and 
more complicated orbits. This finding results from a detailed 
study of these orbits guided by return times into the high 
velocity peak at the sightline I — 6 degrees. The return times 
display several peaks that are associated with the different 
orbit families (and some scatter by stars missing the ob¬ 
servation window) and help with the quantification of the 
contribution of different orbit families. We note that, in re¬ 
sponse to our work, Molloy et al. (20151 have updated their 
paper and now also conclude that higher order orbit families 
(they name 3:1 and 5:2) can make significant contributions 
to the high uios features. 


Another interesting statistic is the ratio of the number 
stars in the high uios region of a certain sightline divided by 
all stars passing through the sightline at arbitrary uios values 
for all stars that have been observed once in the high velocity 
peak. While we find that xi orbits with stable short return 
times show an average ratio of u ~ 0.45, the stars in the 
simulation start with a ratio of ~ 0.22, as a consequence 
of the significant contribution of more complex orbits to the 
velocity peak. 

Within 2 Gyr of the initial detection, this ratio drops 
steadily down to iz ~ 0.12. Due to changes in the bar poten¬ 
tial, length, and pattern speed, as well as orbital diffusion, 
stars diffuse out of these orbit families weakening the high 
uios peak. However, the high uios feature in our simulation 
persists, as the growing bar captures new stars from the 
surrounding disc. It is unclear, how realistically our model 
represents the precise balance of orbit capture and diffusion, 
which is determined by chaotic diffusion, two-body interac¬ 
tions and fluctuations in bar parameters. In addition, obser¬ 
vations of enhanced star formation at the bar tips and on 
bar orbits along the leading bar edges (see above), imply a 
source of young stars on high velocity bar orbits, which is 
not mirrored by our model. Overall, our result provides an 
indication for the recent growth of the Galactic bar. Quan¬ 
titative constraints on the bar dynamics and history would 
demand further exploration of model parameters and also a 
detailed age and distance tomography of the observed high 
Wios stars. 


We note that our bar is a slow bar with 7?bar ~ 0.5i?cR 
in a disc with a high dark matter contribution to the circular 
velocity curve. As discussed in Section 4, observations do 
not constrain the pattern speed and the length of the MW 
bar well and we currently can exclude neither a slow bar 
as ours nor a fast bar with f?bar = 0.7 — 1.0i?cR- As far as 
the dark matter fraction is concerned, our disc agrees with 
recent estimates on the local circular velocity in the MW 
and the corresponding DM contribution. Further work is 
necessary to understand how faster bars and/or lower dark 
matter fractions would alter the processes discussed in this 
paper. 


Debattista et al. (20151 point out that a nuclear stellar 
disc with a radial extent of i? ~ 1 — 1.5 kpc could explain 
high velocity features at Z ~ 10 degrees. Although this idea 
is interesting, they do not currently have a model which can 
explain the velocity distributions at the full APOGEE range 
of bar longitudes, as our model can. The simulations of|Sor-| 


mani et al. (20151 advocate against such a giant X 2 disc for 


the current MW centre, but the extent of such a disc may 
have varied over time. To safely determine whether this idea 
can explain the data, a model, which self-consistently pre¬ 
dicts velocity distributions at longitudes I = (—20) — 20 de¬ 
grees would be needed. A higher number of observed stars 
in this region with reasonable measurements of distances, 
ages and metallicities would certainly help to (in)validate 
our model or other attempts to explain the observed veloc¬ 
ity distributions. As far as the nuclear stellar disc in the MW 
is concerned, its existence so far has only reliably been de¬ 
tected in the central 150pc ( Schonrich et al.||2015 ( and the 
short scalelength inferred for this structure currently seems 
to disfavour a large nuclear disc as proposed by |Debattista| 


et al. (20151. 


To summarise, we have shown that the high viob fea- 














































tures seen towards the bar in APOGEE are likely due to 
preferentially young stars on orbits moving along the major 
axis of the bar. They move close to the plane, which is why 
the feature disappears at high latitudes. The most impor¬ 
tant orbit family contributing to the feature are xi orbits, 
but a variety of more complex orbits also make significant 
contributions. Despite its simplicity, our model is the first 
shown to match the detailed kinematics in the bar region, 
and hence its late stages can serve as a blueprint for further 
and more detailed studies. 
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